
function Complete_Matlab(option)
global p  q  dt  tmax  y ;
y0=[p;q];

dt=1;

tspan = 0:dt:tmax;  % time vector,
options1=odeset('RelTol',1e-6,'Stats','on');
options2=odeset('RelTol',1e-8,'Stats','on');

tic
if option ==1
    disp('ODE 45 Explicit Method Selected');
    [t,y]=ode45(@waves_benchmark,tspan,y0,options1);
else
    disp('ODE 15s Implicit Method Selected');
    [t,y]=ode15s(@waves_benchmark,tspan,y0,options2);
   %[t,y]=ode15s(@upwinding,tspan,y0,options2);
    
end
toc


